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Abstract 

We present analytic results for the polarizabilities, x{Qt^)j iii the charge, 

spin, and current channels for two-dimensional s- and d-wave superconductors 

for large momentum transfer. A collective mode in the charge channel is 

predicted to exist for extremum vectors of the Fermi-surface with energy below 

twice the maximum superconducting gap. Such modes are directly observable 

through inelastic x-ray or electron scattering. Scattering of single particle 

excitations by these collective modes leads to several unusual features in the 

single particle spectrum in the superconducting state which are seen in angle 

resolved photoemission experiments. 
PACS numbers: 74.20.-z, 74.20.Fg, 74.25.Nf 
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I. INTRODUCTION 



It is now well established that the hole doped cuprates superconductors are d-wave 
superconductors with a very strong two-dimensional anisotropy. [0] It is clear that electron- 
electron interactions lead to the attractive pairing, although no consensus has emerged 
about the precise physical nature of such interactions. These interactions must also produce 
a very unusual spectra of particle-hole fluctuations leading to the observed non-Fermi-liquid 
transport properties in the normal metallic state, at least, near the hole-density for maximum 
T 

A phenomenological model |^ from which many of the anomalous transport properties 
can be calculated consistently assumes that the particle-hole fluctuations are scale-invariant, 
as in the critical region of a quantum phase transition. This model further assumes that the 
fluctuations are nearly momentum independent, leading to the prediction that the single- 
particle scattering rate - measurable in angle resolved photoemission(ARPES) or tunneling 
spectra - has the same anomalous temperature (and frequency) dependence as in the momen- 
tum scattering rate measured through the resistivity or the energy scattering rate measured 
by thermal-conductivity. The prediction of an unusual single particle line shape was verified 
by ARPES experiment. ||^ 

Any particle-hole spectra, usual or unusual, must develop a gap related to the super- 
conducting gap below Tc. If the lifetime of quasi-particles in the normal state is primarily 
determined by electron-electron interactions, it must become very long in the superconduct- 
ing state compared to the extrapolation from the normal state. [Q This prediction has been 
also verified. A further prediction followed from the realization that the particle-hole 
spectra must revert to its normal state form at energies well above twice the (maximum) 
superconducting gap. Therefore quasi-particle line shape for energies above about 3A must 
continue to be of the normal state form. The single-particle spectrum in the superconduct- 
ing gap must then be a sharp peak near the BCS energy followed by a broad hump starting 
at about 3A. These predictions have been also qualitatively verified through ARPES 
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spectra. 



Recent more precise ARPES measurements |Ty,|Tl[ have motivated us to examine more 
precisely the change in the particle-hole spectra and the single particle spectra in the super- 
conducting state. The ARPES spectra shows a slower variation with ek of the sharp peak 
than BCS and a more robust hump near 3A than expected. Also, the calculations must be 
performed using the d-wave anisotropy of the gap rather than the s-wave gap used in the 
simple calculation performed earlier 

To obtain the single particle spectra, we must calculate the polarizabilities for the charge, 
current and spin channels as a first step. We expect that the momentum-dependence of the 
polarizability of d-wave superconductors would be qualitatively different from that of s- 
wave superconductors. The result for the polarizability of s-wave superconductors, at large 
momentum for both two and three dimensions (and the associated anomalies in the phonon 



spectra), were presented earlier. |1T2|-[14| We include the results for polarizability of s-wave 
superconductors, to compare with d-wave superconductors. For a spherical Fermi-surface 
in three dimensions, the lowest order polarizability in charge channel has a logarithmic 
singularity at energy 2A for momentum close to the extremum vector spanning the Fermi- 
surface, {2kF, 0, 0). In two dimensions the logarithmic singularity is present near momentum 
(kp, kp) but for momentum close to (2/ci?, 0), the singularity is stronger ~ (cu — 2A)~^/^. We 
give further mathematical details in the Appendix A. 

The singularities in the superconducting state in the polarizabilities at large momentum, 
(while they are absent in the normal state) arise because the singularity in the one particle 
density of states in BCS-superconductors near A are strong enough to lead to singularities 
in the joint particle-hole densities of states at special momentum-vectors connecting points 
on the Fermi-surface. In the d-wave case, these will be vectors connecting points on the 
Fermi-surface with the maximum gap. Taking these points to be along the x and the y-axes 
in two-dimensions, see Fig. 1, we can define two vectors Qo and Qa- 

kl 

Ck.Qa = -vpkx + 
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Qo will often be called the extremum vector. For an isotropic Fermi-surface, or more gener- 
ally with a Fermi-surface of a square lattice we can interchange x and y in Eq.(|I]). 

We find that no nesting is required for the singularities. However nesting, i.e., ^ = 
for points of the Fermi-surface (coincident with the maximum of the gap for the d-wave case 
) leads to stronger particle-hole singularities. 

The most important new result in this paper is that the particle-hole spectra in two- 
dimensions has a singularity at energy below 2A at large momentum transfer in charge 
channel in d-wave superconductors leading to new collective modes. Such collective modes 
can be observed directly by inelastic x-ray or electron scattering. Thus the particle-hole 
spectra does not just develop a gap below 2A as earlier assumed. These results are true 
even when the normal state is a Fermi-liquid. 

The collective modes also strongly influence the single particle spectra measured by 
ARPES. We compute the modified single-particle spectral function through the self energy 
near the maximum gap. We find a sharp peak near A coming from the original BCS form 
and the broad hump below 3A arising from the exchange self energy due to the collective 
mode. 

This paper is organized as follows. We first define the model and show the general formula 
for various response functions in section II. We also discuss the effect of the coherence factors 
in this section. We present analytical(and graphical) results of calculations of polarizabilities 
for both s-wave and d-wave superconductors, in section III and IV, respectively. In section 
IV, we also show single particle spectral function and compare with the experiment. In 
section V, we discuss the predictions, limitation of our results and comparison with existing 
experiments. 

II. MODEL AND DEFINITIONS 

We consider the simplest possible model for the electronic energy 



where /i is the chemical potential, kp/{2m), in a circular Fermi-surface of radius kp shown in 
Fig. 1. Any model with nesting on the Fermi-surface will have more singular polarizability 
than the isotropic surface, Eq. (j^). We shall mention the expected corrections for such a 
case. 

The single particle Green's function in the superconducting state can be written in the 
Gorkov-Nambu notation as [T^ 



G{k,uj) = -2 : , (3) 



where = (e^ + A^). and Ak = A for s-wave superconductor and Ak = A cos (20) for 
d-wave superconductor. 

The lowest order polarizability in charge, spin and current channels at momentum Q 
and energy u are given at T = by 

xt'^'-'^Q^co) = -^J ^r^Tr[r,G{k+Q,u + uy,Gi^,u;% (4) 



(27r)= 

xnQ,^) = -^J ^j^Tr[G{k + Q,u; + u;')G{k,u;')]. (5) 
^<^^ip-\Q^u) = -i—J——{k + q/2U^ (6) 



(27r 

where Tj are Pauli matrices. 

After performing the frequency integral, they are given by 
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charge^r^ n ^ I i ^k+QCk - Ak+Q Ak \ / 



uj + Ei,+Q + El, + ir] io - Ei^+Q - Ek + ir] ^ 

(7) 



° ' 2 ^ \^ ^k+Q^k J \u + Ek+Q + Ek + ir] u - E^+q - E^ + ir] ^ 

(8) 



current 
Xo,ij 



(Q, oc 1 ECk + Q/2),(k + Q/2), (l - !t±Q^±^i±Q^) 
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(9) 



\uj + Ek+Q + Ei, + ir] to - -Ek+Q - E]^ + ir]J 
We are interested in the response at large Q, comparable in magnitude to kp- In such a 
case, in a lattice, the current response always has a part proportional to the charge response. 
PJ27[| This is missing in the translationally invariant model of Eq. (j^). 



Let us discuss the importance of the coherence factors. Note that the coherence factor 
contains — Ak+qAk in the charge channel, while +Ak+qAk in the spin channel. Therefore 
in the case of s-wave superconductors, any interesting features near Fermi-surface can arise 
only from the charge channel. However, d-wave superconductors have a change of sign of 
the gap under 7i/2 rotation: 

A^ = -A<^+,/2. (10) 

This leads to the result that the susceptibility at Q ~ {kp, kp) from antinode(-l-) to 
antinode(-) is negligible in the charge channel, while the susceptibility at Q ~ {'2^kp, 0) 
from antinode(+) to antinode(+) is negligible in the spin channel. Therefore, the interest- 
ing features in charge channel can arise from the excitation of Q ~ (2A;^,0) and that in 
spin channel from the excitation of Q ~ {kp, kp) near uo = 2A. This issue is related to the 
sharp magnetic neutron scattering peak which was observed in the superconducting state 
of Y Ba2Cu^0-j near momentum Q = (vr, vr, tt) ( which in our notation is {kp^ kp) and odd 
with respect to inversion between the two planes per unit cell), and energy uj = AlmeV. 



Tq-ff8l It was pointed out in Ref. [16| that the BCS coherence factor in neutron scattering 



amplitude vanishes unless Ak and Ak+q have opposite signs, because magnetic scattering 



is odd with respect to the time reversal symmetry. This was argued in Refs. |T6|,|T9|-Q to 
show that the peak is a manifestation of the d^^-y-i state. Another suggestion was that the 
superconducting order parameter has s-wave symmetry and opposite signs in the bonding 



and antibonding electron bands formed with a CU2O4 bilayer. pH] A numerical calculation 



shows a collective antiferromagnetic excitation for momentum (tt, tt, tt) with d-wave gap. 
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26| . These calculations concentrated only on the response function in the spin channel. 

We are interested in the polarizability for the momentum vectors connecting the maxima 
of the absolute value of the gap. Near such a maximum gap, we can expand the gap in 
momentum space, e.g., 

Ak+Qo/2 = A(cos^ (p - sin^ 0) 

'^(kp+Kf+ki '^^^ ^ki'- 

Qo connects the (+) antinode and (+) antinode and Qq connects the (+) antinode and 
(-) antinode as shown in Fig. 2. The dispersion relations, e.g., Eq. (|l]), show that one 
direction of the dispersion is more slowly varying than the other; it is linear in kx{ky) but 
quadratic in fcy(fc^). At a general Q, the dispersion is linear in both k^ and ky. A singular 
behavior of polarizability, if it exists, is therefore more likely for momentum near Qo and 
Qa at frequencies near 2A. 

III. S-WAVE SUPERCONDUCTORS 

The long-wavelength collective modes for s-wave superconductors have been discussed 
earlier. [^,^,^ We concentrate on the results for large momenta. The results for the charge 
polarizability in this case have been already presented. [1^ The details of the derivation are 



given in Appendix A. Near the extremum vector, Qo, in two dimensions, we found 

A 3/4 AT 1 

r 5uj <o 

Imxr^m^)- .A3/.^, 1 (12) 

where Suj = uj — 2A(1 + ^^5'^/4), and A^q = with v, the area of unit cell, (from now on, 
set V = 1.) Here, q is defined to be parallel component of (Q — Qo) satisfies << 1, where 
^ is the coherence length of the superconducting gap, ~ hvp/A. 

Another singularity has been also found for momentum Q^ = {kp, kp) at energy 2A in 
agreement with earlier results(see, e.g., [p5|] ). 



5uj <0 

Wr^^(Qa,-)-^ . . (13) 

2 Ef 



6uj > 0, 



where 6uj = u — 2 A. 

For completeness, we also give the results for three dimensions. Singular behavior is also 
found in three dimensions, but it is less singular than that of two dimensions. 

■D charge, \ , Ep 

Rexo (Qo, - — ^ In r— — 

Ef ^J\5uJ\A 

where iVo = The singularity in the charge channel couple to phonons. Eq. ( ]T^ ) has 
been used to understand the anomalies generated in the lattice vibration spectra in the 



superconducting state in some compounds. [|l4-[14[| It is clear that for nested Fermi- surfaces, 
the polarizability become more singular. We found that the nesting along the y-direction 



for Qo gives 1/J\uj — 2A| in the polarizability. 



IV. D-WAVE SUPERCONDUCTORS 

The calculation of the polarizability in d-wave superconductors is more complicated than 
for s-wave superconductors. Details of the calculation are given in Appendix. B. We present 
the results below. 

We again concentrate on the case of large momentum transfers, near Qo and Qa. Refer- 
ring to Fig. 2, we find that near Qo, at Q = {2kp — g^,, qy) for \qy\ < -y/^ and \qx\ < (-^), 
where fci? = 1, the lowest order charge susceptibility is given by 



-BKiX), -Aql - ^ < 5^ < 4^ (15) 

-5/2, 5u: < -Aql - 
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ro, 



A3 



(16) 



fln(X), 6uj<-Aq^- 



2 J^Wv 



16A ' 



where 



X 



A 



B 



( ./ll + (E2g2)/(8A2^ 



1 -5cjE2/(4A3 



V 



(17) 



and (JciJ = 2A — 2Ag^ — uj. K(X) is the Elhptic integral of first kind. The results above 
show that there is a logarithmic singularity aX u = 2A — 2Ag^ — A/S? / E'jp in Rexo'*"^^^ 
and u = 2 A — Ag^ + (Epq^) / [16 A) in Imx'o'"'^^'^ ■ The susceptibility is plotted in Fig. 3. 
The singularity arises just below 2A. The shift from 2A, of order A'^/Ep, comes from the 
curvature of the gap. In the case of s-wave superconductor [|l^], the singular behavior of 
\6u\''^/^ at Q = i2kF,0) remains at = 2 A. For d-wave superconduictors, the momentum 
dependence of the gap does not wipe out the singularity, but reduces it to logarithmic and 
changes the position of the singularity. 

As explained earlier, the singularity in the polarizability in the spin channel can arise if 
the momentum transfer connects opposite signs of the gap. The susceptibility in the spin 
channel near Q^, at Q = (kp — q^, kp — qy) for Ig^,! < and \qy\ < is given by 



Imxr"(Q,^) 



A^(0)A , Ef 



=ln 

Ep^l + ql ^\5uj\A 
5uj <Q 



Af(0)A 



2 



5uj > 0, 



(18) 



where 6lu = 2A - 2A(g - SA^ql/E% - u. 

As we see, there is also a logarithmic singularity in the spin channel. But it is important 
to note that the amplitude of the singularity of the charge channel is bigger than that of the 
spin channel by the 0{Ep/A). 
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A. Susceptibility in the Random Phase Approximation 



The next step is to see how the lowest order polarizabihty is modified due to the inter- 
actions among the quasiparticles. We present the approximate renormalized polarizabihty 
in the present section. In the random phase approximation, 



^ irk '^'^ 



where we assume that the interaction Vq = —V. This has a singular part in charge channel 
with a pole at frequency loq with weight rg and a regular part which is approximately a 
constant for uj > 1.5 A. The position of the pole, uq, in x is given by 

cuQ = 2A - 2Aql - 4AV^f - 8Ae-^ 

= 2 A - 2Aql - 4AV^f - 8^e ^ . (20) 



The weight of the pole is 



8A _^ 

-e Av 



16v^A Jl + (E2g2)/(8A2) v/B+(i.|..g)/A2 



(21) 



Four parameters, NqV , A/Ep, A/V, and qy determine the position and the strength of the 
poles. Let us fix NqV ~ 1 and A/Ep ~ A/V ~ 0.1. At qy = 0, the position of pole is at 
u = 1.49A and the weight is r = 0.13. At qy = 0.2, the position of pole is aX u = 1.62A and 



the weight is r = 0.087. For the reasonable condition \qy\ < ^A/Ep ~ 0.3, the pole is at 
about ujq = 1.5A and the weight is ~ 0.1. 



B. Polarizabihty with a marginal Fermi-liquid vertex 

The anomalous normal state transport properties near the doping for maximum can be 
understood through a phenomenological scale- invariant particle-hole spectrum, (Eq. (p2[), 
see below) which at frequencies large compared to temperature is a constant, independent 
of momentum(except in the hydrodynamic regime) and of frequency up to a cut-off Uc- 
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lmxA4FL{q, 1^) ~ Ar(0); T « z/ < z/„ (22) 

The polarizability can in general be represented as in Fig. 4. If we take for F a phenomeno- 
logical, momentum independent vertex whose imaginary part is ~ zA^(O)^^ for T << v < Uc, 

XMFLiq,i^) = Tiu)[Y,Ak,gGik + q,ij + u)Gik,u)] , (23) 

\k,ui ) 

Eq. ( p2D is then essentially reproduced. Here G is the normal state Green's function and 
A's are appropriate vertices. In the preceding section, we put F ~ 5q^qi5u],u)' in Fig. 4. This 
leads to major difference between Xmfl in the normal state and Xo of Eqs. (|) and (|^). But 
for the superconducting state, the results are qualitatively similar if we make the further 
assumption that F(i/) does not change in the passage to the superconducting state. This 
assumption is almost certainly incorrect. But as long as F(z/) does not develop any sharp 
features for energies below 2A in the superconducting state, the results are unlikely to be 
qualitatively altered. This is because of the sharp drop in xq for energies below 2A which 



we have calculated. Using ([15D in (0), we see that the log singularities near 2 A change 
to log'^{- ■ ■) singularities. Within the approximation of the phenomenology this difference 
is not important. There is one important difference however, x in Fig. 3 and in Eq. (p3D 
is the total polarizability, not the reducible polarizability. So no further manipulation as in 
section (A) is now not permitted. 

C. Single particle spectral function in the superconducting state 

The singularity in the polarizability at large momentum provides important modification 
in the single particle spectra. The singular contribution to the polarizability in the spin 
channel and in the current channel for translationally invariant case is 0{A/Ef) smaller 
than the charge channel. As we already mentioned, in a lattice the polarizability in the 
current channel is proportional (and usually of similar magnitude) to the charge channel. 
We may therefore ignore the spin channel. 

Generalized self energy can be obtained using the diagram shown in Fig. 5. [|l^ 
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where k — (k, a;). Here we neglect the other self energy contribution shown in Fig. 5(c) ; 

It is clear that this is negligible because it is of the order of 0(ek/-Ek) and we are near the 
Fermi-surface. 

We write x(Q,a;) as the sum of a singular part and a regular part. 



x(Q,a;)^ ^^J^ /(Q)+X.e.. (26) 



We will define S'''"^ as the contribution of the first part and S*"*^^ as the second one with 

E{cu, k) = J:'"'3{u, k) + S"^5(u;, k). (27) 

For a d-wave superconductors f{Q) gives the momentum limits shown as above, Q = 
{2kF — qx,Qy) where \qx\ < A/Ep and \qy\ < ^ A/Ep setting kp ~ 1. We assume that 
cuq — cuq, and tq = r in this range. 

The self energy, E*'"^(a;,k) in the charge channel is 

E--«(u;,k) =i I ^0R{VT^)Go{u; + u;^,k + Q){VT^)x{u;,,Q) 

. f duJi(PQ {uj + uji)l + ek+QT3 + Ak+QTi ^ 2uJor 
J (27r)2 {u + uiy - E^_^Q + ir] 'cuf-cu^ + irj 

After performing the frequency intergral, we found that the singular part of the self energy 

is given by 

d'^Q 1,, 1 

V f 

1 



a; - Wo - -Ek+Q + if] 
E.,(w,k) = V'r I ^{^^){ 



(29) 



■ UJQ - Ek+Q + IT] 
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IT] 

), (30) 



1 



+ %^ )• (31) 

After expanding -Ek+Q near maximum gap(+), we obtained the self energy for ~ 
—kp + k'^ where \k'J\ < \qx\ and \ky\ < \qy\, 

^ ^ , , , V'Nor ^ . 2A , 4A2 
Ke2^/(a;, kj = ^[arcsm U[—uj + 00 ——) 



4^2 JEp(uj* - uj) E_ 



F 



■ 2A , 4A2 

arcsm , + oj 

IEf{uj* + u:) 



T V r Yll^^\ ■ 2A 

imLr cij, k = ^ — arcsmh — U[uj — uj 



2A 

+arcsinh^ e{~uj - uj*)\ (32) 



R V r ^'^0^ • 2A , 4A^ 
KeL^-^^ciJ, k) = ^ [arcsm t^(— u; + 

4V2 \ Ep[uj* -uj) Ef' 



2A ^, , 4A\, 

+ arcsm , U[uj + u ^E^n 

'Ef{uj*+uj) Ef 



^ T ( ■ 2A 

imL^-^(C(j, kj = ^ arcsmh . =0[uj — uj ) 

4V2 JEf(-oj* + oj) 



2A 

+arcsinh 6{—uj — uj*)], (33) 

'Ef{-uj* - uj) 



where tu* = cjq + A — 4AA;,^ + Ep/A{ky — 2k'^)'^. S^g is basically the same as Cfc/A^S^^, 
expanding the ek+q = -Ck + 0(q) with -Ek+Q = -^k- 

We now discuss Xreg- This depends on whether one uses the form Eqs. (|15D and (|16D 
for Xo or the marginal form. This difference is important for large frequencies, because the 
single particle spectra for the former case reverts to the Fermi liquid form, which disagrees 
with experiments. So we consider the marginal fermi-liquid form. An important issue in 
that case is the form of Xreg for low frequencies in the d-wave superconductor. In a d-wave 
superconductor, the regular part of XregiQ, J^) at low energies is dominated by the particle- 
hole excitations with Q connecting the nodes of the gap function Ak. These determine the 
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lifetime of low energy quasiparticles and the low temperature thermodynamics. But they 
are not involved in the scattering of quasiparticles near the maximum of the gap due to 
restrictions of momentum conservation. For the relaxation of such quasiparticles, we may 
therefore assume a marginal Fermi- liquid form for energies above 2 A: 

with Treg + r = 1. This produces 

ImE^(^^) oc T\uj\e{uj - 3A) + \uj\e{-uj - 3A). (35) 

Fig. 6 shows the Re and Im parts of self-energy for excitations near the maximum gap 
including both the singular and the regular parts of the fluctuations. 
The one particle Green's function can be obtained using 



^ a; — ek + — E/ — Ak — 

y — Ak — Et-j a; + Ck — E^-g J 



(36) 



The spectral function is give by: 
A{u, k) = — ^ ^ ImCii 

TT 

^ -gg^M j^ 0^ - ek + - Ej 

TT (a; - ek + E^g - E/) (a; ek - E^-g - E/) - (Ak E^J^ _^ • 

The spectral function is plotted in Fig. 7 for k near the maximum point of the gap and for 

the representative values, V ~ Ep = lOA. For < a; < 2A where ImE/ = ImE^-^ = 0, it is 

a delta function near a; ~ A and has a broad hump at a; ~ 2.5A. As we change k, we find 

the peak near A shifts slower than the change of ek, itself. The position of the broad hump 

shifts by Aky which can be up to order of A'^/Ep- The broad peak would become smooth if 

the contributions from the other excitations, e.g., node-node excitations are included. The 

excitations from node to node contribute to finite life time of quasiparticle for a; > A which 

are not included here. However, this would not change the qualitative behavior of the broad 

peak. 
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V. DISCUSSION 



We have found a new collective mode for large momentum at energy just below 2 A, which 
can be detected by inelastic x-ray scattering or electron scattering. Since neutrons have a 
negligible cross-section for scattering with electronic density fluctuations, this mode cannot 
be detected by neutron scattering. The calculations, done using an isotropic Fermi-surface 
give a lower limit to the spectral weight of the collective modes. Realistic band-structures 
show some nesting in the region of the maximum gap due to which stronger effects are to 
be expected in experiments. 

We have also calculated the single particle spectral function which is strongly affected 
for states near the maximum of the gap by the collective mode. This quantity has been 



measured by ARPES |TT|. ARPES on Bi2212 near (tt, 0) in the superconducting state 
shows two features. There exists a sharp peak around AOmeV and a broad hump at about 
lOOmeV. Our results also show two peaks in the spectral function, a sharp peak near A and 
a broad peak near 2.5A. When we change the momentum by k in a direction perpendicular 
to the Fermi-surface from the point of gap-maximum, the sharp peak and the broad peak 
shift as Epk'^ where k'^ can vary up to A/Ep (See Fig. 7 (d), and Fig. 2 (c) in Ref. [0). 
For a change of momentum transverse to the Fermi-surface, this variation is cancelled to 
leading order by a similar decrease in the value of the gap. (See Fig. 7 (b), (c), and Fig. 2 
(b) in Ref. [jll|) These features are generally consistent with the experiments, although the 
persistence of the sharp feature near A is over a wider region than calculated by us. We 
think this is probably due to the nesting and Van-Hove dispersions of the band-structure 
near the points of gap-maximum. 

It is also apparent in the experiments that when the sharp peak disappears, the broad 
peak start to move rapidly as in the normal state. This is also consistent with our expectation 
that the behavior of spectral function must be similar to that of normal state for Ck away 
from the chemical potential by much more than A. 

We, finally, mention the limitation of our results. Due to the use of an isotropic Fermi- 
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surface, the range of momentum where our results are vahd is rather hmited. Another 
hmitation of our results is that the one particle Green function G is not calculated self- 
consistently. We start with G of the BCS form. The singular self energy corrections split 
the BCS peak in the spectral function into two peaks. But we do not go back and calculate 
the particle-hole fluctuations and the self-energy with the renormalized G. We expect that 
using renormalized particle G, additional features at higher energy( near 5A) appears and 
the features already calculated are sharpened. The magnitude of these higher order effects 
is expected to be of order N{Q)V/\/Ef- We have not included the vertex corrections in a 
systematic fashion either but have contented ourselves with a phenomenological requirement 
that the vertex be such that the spectral function in the superconducting state reverts to 
that of the normal state at large enough energy. Another limitation is the assumption of 
the sharp change as a function of v in Eq. (|3^). We expect a rapidly varying but smooth 
function. This would tend to broaden the onset of the hump in the spectral functions in 
Fig. 7. 



APPENDIX A: 



We show how to evaluate the integrals to obtain the results quoted in Ref. ||12[ and 
section III of this paper. After performing the frequency integral in Eq. one must 

evaluate 

^a^ EkEk+Q - e^ek+Q + ^ 1 ^ 1 ^ ^^^^ 

Ek-Ek+Q -Ek + -Ek+Q -u + ir] Ek + Ek+q + u - ir] 

The coherence factor does not play an important role for large Q so that we simply put it 
0(1). Both i?k and -Ek+q are bounded from below by A. Now expand E'k and -Ek+q near 
the chemical potential in cartesian coordinates and (without loss of generality) take Q along 
X-axis. First consider Qq = 2kpx, 

2 2 

ek = -vpp. + + ^ + 0{pI) (A2) 
2m 2m 
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where px = {kx — kpx), Vy = {ky — kpy) and Pz = {kz — kpz), p's are -C fc^'s, where 
kp = {kpx, kpy, kpz)- (For this Qo, kpy = kpz = 0) 



2 2 

'^^Qo=vfP. + ^ + ^ + 0{pI) (A3) 



Now i?k and -Ek+Q can be expand as follows. 



^^ = ^ + f 

Sk+Q,, = A + 1^ (A4) 



The integral of Eq. ([Al| ) becomes 

1 







^(2A - + ^[t;2p2 + (^2/3^ + p2/2m)2] ^^^^ 



One can see immediately by power counting, for three dimensions, we have log\uj — 2A| 
singularity as p appproachs (means k — > ki?). For two dimensions, we have [u — 2A]~^/^ 
singularity. 

For general Q, again choose Q parallel to x-axis without loss of any generality and define 
Q = 2QxX = 2kFxX (See Fig. A2). We again expand in terms of the p's defined above. 

Qx , kpy , kpz pI pI pI 

= Px + —Py + Pz + 7^ + 77^ + 7^ 

m m m 2m 2m 2m 

Qx kpy kpz pI pI pI . 

Ck+Q - — Px H Py H Pz + 7^ h 7^ h — (^Abj 

m m m 2m 2m 2m 



Unlike Eqs. ( [A2D and ([A3|) , linear terms in Py and Pz are present in Eq. ( [A6|) . This 



changes the power counting. When we do the same expansion as in Eq. (|A4|) and do the 
integral, we see by just power counting that in the lowest order of polarizability the log is 
cut off by the kpy/2m and k'jp^/2m terms in three dimensions. (In fact, it has log singularity 
in two dimensions with coefficient depending on Q and changing to the stronger singularity 
mentioned above at Qo = 2kFx) 

To summarize, singularity at = 2A for = 3 is found for only Qo = 2kFX (or << ^^^) 
where the leading terms of px is linear and Py, pz are quadratic in Eq. 
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APPENDIX B: 



We show here how to evaluate the integrals to get the results quoted in section IV. We 
start with Eq. (|^ in section II. Expand ek near the Fermi-surface and Ak near the maximum 
gap(See Eq. and ([TT|) in section II). For k = {—kp + Px,Py) 

Ek ^ A - 2Apl + ^{2pl - 2p,pl + pI/2), (B1) 

where Ep = kp/{2'm) and Px{y) = Px{y)/kF (we set kp = 1). Note that it contains 2Ap^ term 
due to the anisotropy of the gap. We can also find -Ek+Q,, ii^ar Qo; {k^, ky) = {kp — Qx + 
Px,Py + Qy)- The expression for E'k + -E'k+Qo — is written as follows. 

f{Px,Py) = Ek + Ek+Q - UJ 

= ^{^pI + ^P^QyPy + + + ' " 4ApJ + 2A - 2AgJ - uo. (B2) 

The susceptibility is obtained by integrating the 1/ f{px, Py) with respect to px and Py. After 
integrating over p^ {—A/Ep < p^ < A/ Ep) we have 

A 1 



J ^ El + ^^^^ _ + ^ _ ^ 



X arctan , ^^^^ ^ ^ = (B3) 



The leading term of this function is 

TT A 



(B4) 



After performing the intergration over py, we have the Elliptical function. 

K{k) = f '^'^ (B5) 







where 



(B6) 



When two poles of Elliptical function coincide( k = 1), we have a log singularity, which 
happens here at 6lj = AA^/Ep. 

1 f 1 



ReXo(Qo,^) = -^/y^c^Mp. 



/8 + (E2g2)/A2) VI'^^-4AV^^| 

where iVo = l/iiTrEp). 

The singular part of self-energy for spin-spin correlation function is also presented here 
which are similar to those of the charge channel apart from signs. 

KeL/ = ^ arcsm , U[—uj + 6 A -— ) 

2V2 JEf{3A - u) Ef 



2A ^, , 4A^, 

— arcsm , Oioj + 3A ;=— ) 

lEpi'iA + uj) Ef 



ImS/ = ^ — arcsinh^^ 9{uj — 3A) 

4V2 JEf{-^A + uj) 



2A 

+arcsinh^ -.e{-uj - 3A)] (B8) 

lEpi-'iA-uj) 



V^N^r^ 2A ^ 4A 

ReS^j = — : — T^larcsin — , 0(—uj + 3A 



2 



2V2 ./^^(3A - uj) 



2A 4A^ 
+arcsin 9{uj + 3A )] 

/e^(3A + c<;) Ef 



ImS.. ^ -^[arcsinh , Oj. - 3A) 

4^2 JEf{-3A + u;) 



2A 

+arcsinh^ -.Oj-uj - 3A)] (B9) 

'Ef{-3A - uj) 



Note that the weight, r, here is smaller than that in Eqs. (^) and (|33D by order of A/Ef- 
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FIGURES 

FIG. 1. The isotropic Fermi-surface with radius kp illustrating the vectors, Qo = {2kF,0) and 
Qa = (kp, kp), used in the text. 

FIG. 2. The isotropic Fermi-surface with d-wave gap symmetry illustrating vectors defined in 
the text. The and (-) denote the the signs of the gap. Qo and Qa are momentum transfer 
between two maximum gaps with same and different signs, respectively. The range of qx and qy 
where our results for the polarizabilities are valid are shown as shaded ellipse. 

FIG. 3. (a) The imaginary part of the lowest order polarizability, ImxQ'"''^^(Qo, w). (b) The 
real part of the lowest order polarizability, RexQ'*"'^^'^(Qo, to). The parameters are the same as those 
in the text; NqV ~ 1 and A/Ep ~ A/F ~ 0.1. 

FIG. 4. The diagrammatic representation of the polarizabihty with vertices A and 

FIG. 5. (a) The renormalized interaction after RPA. (b) The lowest order self-energy; exchange 
contribution from the Coulomb interaction (c) The lowest order self-energy; direct contribution 
from the Coulomb interaction. 

FIG. 6. The lowest order self energy, (a) ReS/(cj,k), (b) ImS/(u;,k), (c) ReS^-i (a;, k), and (d) 
ImSriiio, k), where k is near the maximum gap. The parameters are the same as those in the text; 
A^o^ ~ 1, A/Ep ~ A/V ~ 0.1, and r ~ 0.1 

FIG. 7. The spectral function ^(a;,k) with (a) k = {-kp,0), (b) k = {-kp,0.1kp), 
(c)k = {—kp,0.2kp), and (d) k = (— 1.05A;f, 0). The parameters are the same as those in the 
text; A^o^ ~ 1, ^/Ep ~ A/V ~ 0.1, r ~ 0.1, and rreg ^ 0.9 
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